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FOREWORD 

Work described in this report was accomplished by the Flight Mechanics 

Division of the Air Force Flight Dynamics Laboratory and the Digital 

programming section, 4950th Tecc Wing under Project 1431, "Flight Path 
Analysis," Task 143109, "irajectory and Motion Analysis of T light 
Vehicles." The forr-jiation and interim documentation were completed 
by Major Urban D. Lynch. Programming was accomplished by Mr. Fay 0. 
Young of the Digital Programming Section (ADDP), Computer Science Center 
4950th Test Wing. 


This report was prepared by Mr. John J. Dueweke of the High Speed 
Aero Performance Branch (FXG), and Mr. Fay 0. Young, and combines the 
applicable portions of FDL-TDR-64-1, Part I, Volume 1, with vhe interim 
documentation. The overall report is divided into four parts: 

Part I. Capabilities of the Takeoff and Landing Analysis Computer 
Program 

Part II. Problem Formulation 
Part III. User's Manual 
Part IV. Programmer's Manual 

This report was submitted by the authors in March 1972. 
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ABSTRACT 

A well-defined integration of the various aspects of the aircraft 
takeoff a’nd landing problem is presented in the form of a generalized 
computer program. Total aircraft system performance is evaluated during 
the glide slope, flare, landing roll and takeoff. 

The flight dynamics of a generalized, rigid-body, aerospace vehicle 
are formulated in six degrees of freedom. A flat, nonrotating Earth is 
assumed. The independent equations of motion of up to five oleo-type 
landing gears are also formulated. 

A control management formulation is developed to automatically adjust 
control variables to correct errors in the "ehicle's dynamic state. 
Stability in the small is used to maintain stability in the large. 

The equations of motion are integrated using a generalized variable- 
step Runge-Kutta technique. 

The formulation is programmed for the CDC CYBER 70 and 6000 series 
computers using the SCOPE 3.4 operating system. The entire program 
is written in Fortran Extended 
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SECTION I 
INTRODUCTION 

In the design of an aircraft, the engineer is confronted with the 
problem of takeoff and landing and the design of aircraft systems and 
techniques to perform this function. The final evaluation of these 
systems lies in the answer to the question: How does the aircraft and 
its systems perform as a unit? The Takeoff and Landing Analysis (TOLA) 
Computer Program is the result of an attempt to generalize the aircraft, 
the main aircraft control systems, and the landing-takeoff situation 
into a single comprehensive calculation to answer this question. 

The TOLA simulation answers the above question in the form of a 
well-defined integration of the various aspects of takeoff and landing. 
In the equations of motion the assumption is made that the main aircraft 
frame is rigid; however, the dynamic effects of up to five independent 
landing gears are included in the equations. The position and velocity 
of each strut and secondary piston are obtained by numerical integration 
subject to position constraints (for example, the main strut must move 
within the limits of the fully extended position and strut bottoming 
position). The same form of solution applies to the aircraft itself. 
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SECTION II 

COMPUTER AND SYSTEM REQUIREMENTS 

The Takeoff and Landing Analysis Computer Program (TOLA) has been 
written for use with the CDC CYBER 70 and 6000 series computers using the 
SCOPE 3.4 operating system. The entire TOLA is written in Fortran 
Extended. 

1. CYBER 70 or 6000 series computer 

a. A CYBER 70 or 6000 series computer with 32K (decimal), or 
larger core. 

b. Six CDC tape transports. 

c. Control Data Card Reader. 

d. Control Data Printer. 

e. 12 inch CALCOMP off-line plotter. 

2. The CDC tape transports may be replaced by other equipment which 
will simulate magnetic tapes such as disk storage, except for one tape 
unit that may be used to generate a plot tape for the CALCOMP off-line 
plotter. 
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SECTION III 
PROGRAMMING CONCEPTS 

1. THE USE OF COMMON 

Whenever possible, a variable is placed in the FORTRAN "COMMON" Area. 
There are several reasons for this: 

a. The communications between subroutines is simplified. 

b. The structure of the directory is simplified. Since the number 
of variables ir, COMMON is quite large, all COMMON cards are not 
placed in each assembly/compilation. Instead, required "duirmy" 
cards are placed in each deck of source cards. This has in a 
small manner reduced the number of COMMON cards in each deck. 

2. TABLES AND TABLE USAGE 

One of the usual required modifications of any program is the change 
of table sizes. With this in mind, a COMMON block of locations has been 
set aside and the required number of ceils for each tabie is specified 
with data (see TABRE for data preparation). This requires no reassembly 
or recompilation unless the total number of cells required exceeds tne 
COMMON block of 800 cells. 

3. SYMBOLIC INPUT 

Although the FORTRAN system itself has a system of input routines, 
the program does the actual translation of the cards using special coded 
routines. Input data may be read using a system of symbols which is 
designed to give engineering meaning to the analyst. The symbols are 
referenced to actual locations by the use of COhWON and subscripts. 
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4. TRAJECTORY PRINTING METHOD 

The printing of a trajectory may be divided into four categories. 

a. Initial Printing - The printing of specific values at the 
first stage and at each subsequent major stage. 

b. Code Printing - The printing of codes which will identify 

the variables which are to be obtained in the coming time history print. 

c. Time History Printing - The printing of values specified at 
the requested points of the trajectory. 

d. Diagnostic Error Printing - The printing of errors detected 
by the program. 

All input data involved for a case is printed on the output page 
preceding the computation of the first stage printout. Also, data read 
in at stage times will be printed out between the stages of the trajectory 
output. 

Initial print is designed to print certain values which will be 
constant during the trajectory and serves as a reminder of what valuer 
have been used for these constants. 

Code printing is performed once per major stage to identify the time 
history. 

The time history print is designed to print in a minimum space. 

That is, if a certain variable is not desired as output, it is not printed 
and other desired variables are moved in the print format accordingly. 
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The entire printing is controlled to print on a page 11 x 14 inches 
and will print a maximum of 51 lines per page. Page ejection and lines 
control are provided by the subroutines DEF and LINES. 


TAPE USAGE 



Tapes 

Equipment 

Usage 

Tape 5 

Disk or Tape 

Data Input 

Tape 6 

Disk or Tape 

Printed Output 

Tape 13 

Disk or Tape 

Data saved to be used to 



generate a plot tape by 



Plot program (PLTSDF) 

Tape 16 

Disk or Tape I 1 

| Used by the symbolic 

Tape 31 

Disk or Tape ( , 

| input routine to save 


1 

i input data 

Tape 7 

Tape 

CALC0MP Plot tape 


The above describes the tape usage other than for the FORTRAN system. 
All modification of tapes required may be made with control cards placed 
in front of the program before submitting to the computer. 

6. STRUCTURE OF PROGRAM 

Due to core storage limitations (32K), it was necessary to use the 
Overlay feature. The following is the structure of the program: 

a. OVERLAY (0,0) (FOR. EXT ) 


TOLA 

STFL 

TRNPOS 

FXE 

STFLD 

HI HO 

INUPD 

STOVAR 

TLU 

LN'JPD 

ARRAY 

TFFS1 
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INPUZ LINES 


INTEG 

ASIN 


UPDAT 

ACOS 

VPCS1 

MININ 

A TAN 2 

SACS! 

LGDET 

ERROR 

aeroi 

STGTSI 

EXERR 

AOUAD 

STGTST 

NDTLLi 

OPT1 

TDATA 

ATMS 

LGKAR1 


INVR3 

LGEA3C 

ASRCH 

MULT31 

SDFI.GP 

DEF 




OVERLAY 

(1,0) (FOR EXT) 

TOLAN1 

TSRCH 

TFFS2 

DSERCH 

READ 

PACKRR 

DIPLAC 

RITE 

TAB RE 

DIRODA 

READA 

DIR]DA 

STORE 

DIR2DA 


WRCARD 
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c. OVERLAY (2, 0) (FOR EXT) 


TOLAN2 

ENGREV 

AUTS 

ENGFL 

FLARE1 

CENCL 

AUTPR1 

CTENGL 

THAUTS 

TFFS8 

TFFS9 



7. PROGRAM ORGANIZATION 

The TOLA Computer Program is written in FORTRAN Extended. The 
program is segmented and takes advantage of the FORTRAN overlay features. 

This section attempts to describe the overall organization of the 
program from the viewpoints of control cards, tape usage, deck set-up, 
and organization. 

The program is broken up into three overlays as follows: 

a. Overlay (0, 0). Contains all system routines, main, executive, 
integration, computation of the equations of motion, and printing. 
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b. Overlay (1, 0). Set up tables, input routines, and 
input directory. 

c. Overlay (2, 0). All routines of the Autopilot. 

Plotting tapes are generated by a separate program for plotting on 
the Cal Comp plotter. 

a. Storage Reference. All variables requiring arrays have been 
arranged in the standard FORTRAN convention; for example, an array A^ is 
stored in increasing storage locations for increasing i. Matrices are 
stored columnwise. 

b. Integers . All integers are assumed to be in a 60-bit word, 
right-justified. 

c. C^MM^N. In order to decrease the length and time required in 
calling sequences, liberal use of labeled C0MM0N has been made. For the 
actual variable and their arrangement in C0MM0N, the user is referred 
to the program listing. 

d. Variable Names. Because any variable may be referred to by 
FORTRAN, all integer variable names begin with the leading letters I, J, 
K, L, M, or N. This does not mean that all noninteger variable names 
begin with letters other than I, J, K, L, M, or N. They may, in some 
subprogram, be declared integer or real. 

8. DECK SETUP 

a. Running the TOLA Program requires a particular deck setup. The 
deck structure is presented as a guide only in determining this se^up. 
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b. CONTROL CARDS (CDC CYBER 70 or 6000 Series, SCOPE 3.4). All 
control cards are left justified in column 1. The end of record is a 
7, 8, 9 punched in column 1 and on end of job card is a 6, 7, 8, 9 punch 
in column 1. In the control card examples an e.id of record and an end of 
job will be used in place of the cards. 

(1) The following control cards will execute the TOLA Computer 
Program from an UPDATE tape and not print a listing of TOLA. 

Job Card 

LABEL, 0LDPL, R, L = T0LACP, VSN = tape No. RING 0UT 
UPDATE, F. 

FTN, I = COMPILE, L = 0. 

RETURN, 0LDPL. 

LDSET, PRESET = ZERO. 

L0AD, LG0. 

N9G0. 

TOLA. 

End of Record 

Changes to T9LACP if any in UPDATE format 
End of Record 
Data Cards 
End of job. 

If a listing Is desired, omit the parameter L = 0 on the FTN control 
card. 
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(2) The following control cards will generate a new UPDATE tape, an 
absolute file on tape, and list the T0LA Program. 

Job Card 

LABEL, 0LDPL, R, L *» T0LACP, VSN = tape No. RING 0UT 
LABEL, NEWPL, W, L = T0LACP, VSN = tape No. RING IN 
UPDATE, N, F. 

FTN, I = COMPILE. 

RETURN, 0LDPL. 

UNLOAD, NEWPL, 

LABEL, TOLA, W, L * T0LACPABS, VSN = tape No. RING IN 
LDSET, PRESET = ZERO. 

LOAD, LGO. 

NOGO. 

TOLA. 

End of record 

Changes to TOLA if any in UPDATE format 
End of Record 
Data cards 
End of job 

(3) The following control cards will execute from an absolute file 
on tape, TOLA is an absolute file on a tape. 

Job card 

LABEL, TOLA, R, L = TOLACPABS, VSN = tape No. RING OUT 
TOLA. 

End of record 
Data cards 
End of job 

10 
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(4) If it is desired to save data on tape for CALC9MP plotting, 
include the following control card with the LABEL cards: 

LABEL, TAPE 13, W, L = T0LADATA, VSN = tape No. RING IN 

(5) The following control cards will generate a new updated program 
on permanent file (PF), generate on absolute file of T9LA on PF, and 
execute T0LA. 

Job Card 

LABEL, 9LDPL, R, L = T0LACP, VSN = tape No. RING OUT 
REQUEST, NEWPL,*PF. 

UPDATE, N, F. 

RETURN, 9LDPL. 

CATALOG, NEWPL, T0LACP, RP = 999, CY = 1, ID = Prob No. 

RETURN, NEWPL. 

F TN, I = COMPILE, L = 0 
REQUEST, T0LAP, *PF. 

LDSET, PRESET = ZERO. 

LOAD, LGO. 

NOGO. 

TOLA. 

CATALOG, TOLAP, TOLACP, RP = 999, CY = 2, ID = Prob No. 

End of record 

Changes to TOLA if any in UPDATE format 
End of record 
Data cards 
End of job 
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(6) To execute from a permanent file, use the following control 
cards: 

Job Card 

ATTACH, TOLA, T0LACP, CY = 2, ID = XXXXXX. 

TOLA. 

End of record 
DATA cards 
End of job 

(7) To execute from a permanent file, and generate a plot tape 

for the CALC0MP plotter on TAPE 7, use the following control cards. The 
plot tape generating program (PLTSDF) is located on PF: PLTSDF, CY = 1. 

Job Card 

ATTACH, TOLA, TOLACP, CY = 2, ID = XXXXXX. 

TOLA. 

RETURN,, TOLA 

REQUEST, TAPE 7, MT, HI, N, VSN = Tape No. Ring IN 
ATTACH, PLTSDF, PLTSDF, CY = 1, ID = XXXXXX. 

PLTSDF. 

End of record 
Data cards for TOLA 
End of record 
Data cards for PLTSDF 
End of job 


12 


AFFDL-TR-71-155 
PART IV 


9. DATA FORMAT 

Card Format - The program input routine (READ) expects the following 
format. 

. _ I II a ■ • t 

73 - 80 
VIII 

Card Field I - Contains the symbolic name of the variable which data 


Card-Columns - 

1 - 6 

7 

8-10 

11 1 12 - 66 

67 - 72 

Field 

I 

II 

III 

IV 1 V 

VI 


contained in Field V begins loading. Example: Card Column 


1 

IGAM7D 

SIG7D 


12 

-1.23 

90. 
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Card Field II Not Used 

Card Field III - Contains the words DEC, 0CT, BCD, TRA, INT, or is 
blank depending on the type of data to be loaded. The word OCT indicates 
that the data is to be interpreted as octal numbers. The word BCD 
specifies that N binary coded decimal words (N punched in column 12) 
beginning in column 13 are to be loaded. The word TRA denotes to the 
input routine that all data has been read and to return control to the 
calling program. The word DEC and blank are equivalent and specify 
that data loaded is decimal data. 


0CT Example 
Card Column 

1 

8 


12 


NSMAIN 

OCT 


17 

BCD Example 

Card Column 

1 

8 


12 


REM 

BCD 


3SDF2-GEAR-M0D 


The 3 in Column 12 specifies 3 words where each word is considered to 
be 6 characters including blanks. The largest number of 6-character 
words that can be loaded from one card is 9. The analysts should be 
very careful to see that the BCD information does not get punched into 
Field VI. This will cause an input error. 

DEC Example 

Card Column 1,8 12 

VTAB01 DEC 2,0.,1.67,20000.,1.67 

Note that the first character in Column 12 is an integer; the input 

routine will load only one integer per DEC card, and that has to be the 

first number punched in Field V. 

VTABOI DEC 2.,0.,1.67,20000.,1.67 
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If the above card is punched, the two will now be loaded into the machine 
as a binary floating point number. The other numbers will be loaded 
the same, with the decimal point assumed right-justified. 

If anything other than 0CT, BCD, INT, TRA, or blank appears in 
Field II then the word DEC is assumed. 

INT Example 


Card Column 

1 

8 

12 


IP 

INT 

1 


IP 


1, 


IP 

INT 

1,U,1 


When the word INT is useci it is assumed that all numbers on the card 
will be loaded as integers. If only one integer is punched per card the 
INT may be punched or omitted. 

Card Field IV-Not Used 
Card Meld V 

The actual input data to the program is punched in the Field V. 

DEC, INT and 0CT numbers must always be left-adjusted; that is, it 
must always start in column 12 on the input card. All numbers are 
separated by a "comma" and the field terminates with the first blank. 

BCD information begins in Column 13 and the maximum number of 6-character 
words per card is nine. Note that since Field V ends with the first 
blank, the user may punch any comments in the remainder of the field. 

Card Field VI 

This field specifies the initial subscript of the data in Field V. 

If this field is blank, an initial subscript of 1 is Implied. The 
subscript may appear anywhere within the field. 
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Example 
Card Column 


1 

12 

PZER0 

30470.4,41538.24,41538.24 

PZER0 

42538.24,42538.24 


(or blank) 


In the example above, the number 30470.4 is loaded into the first cell 
of the array PZER0. On the second card, 42538.24 is loaded into the fourth 
cell of the array. The one and four punched in Field VI indicate the 
subscript for the array P2.ER0. 


Card Field VII 

Not used as far as the input routine is concerned. This may be used 
as a sequence number for the card. 


10. TABLE FORMAT 

The various types of tables used by the program may be classed as 
follows: 


One Dimensional Tables 

Example 1: NTIRE5 

Card Column 1 

NSTRUT 

NTIRES 


n. • f (i), i = I, 2,..., NSTRUT 

12 

5 

14.,5.,6..6.,6. 


INSTRUT= Fixed point number which is the number of struts on the aircraft. 

For example, the number of tires on strut 2 is 6} i.e., - f(2) = 6. 

i = independent variable values 

n.j = Corresponding dependent variable values 

Example 2: Aerodynamic Data 

Card Column 1 12 

INDA01 1 

ATAB01 .0065,.00748 
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INDAOI 7*1 designates that there are data in ATABOl. The first data point 
is for full ground effect; the second data point is for no ground effect 
in all aerodynamic tables. 


Two Dimensional Table 


Example: 
Card Column 


VTAB01 

1 

VTAB01 


12 


x c C = f W 


N,M 1* X C6 1 ,M z* X CG 2 , ‘** ,M N ,X CG n 

N = Fixed point number equal to 2 times the number of independent variables. 
For a 20-point table, N would equal 40. The total number of machine 
cells required for this table is 41. 

M.j = Independent variable values 

X rr = Corresponding dependent variable values 


N - Dimensional Table 


Example: 
Card Column 


T 

1 


F(N,K ) 
12 N 


IT10W 

NN 

IT10X 

NMN 

TTAB10 

N 1' N 2» N 3 , ' ,,,N NN 

TTAB10 

M N1 ,M N2 ,M N3 , *‘ ,,M NNMN 

TTAB10 

T NrMNl ,T N2 ,M Nl ’•••» T NNN ,M N1 

TTAB10 

T N1 ,M N2 ,T N2 ,M N2. T NNN* M N2 

• • • 

TTAB10 

T N1 ,m nnmn ,T N2 ,m nnmn•••• ,t nnn 


A 


NN and NMN are fixed point numbers of independent variables. , n^,. 




T n , M^^ are values of independent variables. The table subscripts 
NN 

would apply to the N-dlmensional table as well as the two dimensional. The 
total number of machine cells required for an N-dimensional table equal 
NN * NMN + NN + NMN. 
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Examples: 

C = F(X, Y) 

Machine cells required 
C = F(X, Y, Z) 

Machine cells required 


NX = 2 = points for X 
NY = 2 - points for Y 
2x2 + 2 + 2 = 8 cel 1s 


NX = 20 = points for X 
NY = 10 = points for Y 
NZ = 15 - points for Z 
20 X 10 X 15 + 20 + 10 + 15 = 3045 cells 
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SECTION IV 

FORTRAN EXTENDED OVERLAY (0,0) 

1. TOLA - MAIN PROGRAM 

a. Purpose - Initializes parameters in the read routine through 

COMMON, initializes table cata, and parameters in storage routines, and 
calls EXE. 

b. Usage - Calls the executive routine (EXE) for each case. 


2. EXE - EXECUTIVE ROUTINE 

a. Purpose - To zero all variables that may be read in as input, 
initialize subprograms, and set nominal values. Read input, do all the 
proper initialization, set up tables dimensions, check for staging, 
printing, etc. 

b. Usage - The executive routine is the controlling program. When 
a case is completed, a return is made to TOLA. 


3. STGTSI, STGTST - STAGE TESTING ROUTINES 

a. Purpose - To test the possiblity of staging on any of up to four 
increasing and four decreasing variables. 

b. Method - Given the four increasing variable BCD names (in arrcy 
AINCRS) and the four decreasing variable BCD names (in array DECRES) the 
routine first searches the directory for their location (routine STGTSI). 
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FLOW DIAGRAM-MAIN F ?'OG RAM (TOLA) 
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FLOW 01A6RAM - EXECUTIVE ROUTINE (EXE) 




INOSTG 


#0 


TIMEA ■ T + DELTS- I- E- 6 
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FLOW DIAGRAM-EXECUTIVE ROUTINE (EXE) 



i 
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FLOW DIAGRAM - EXECUTIVE ROUTINE (EXE.) 



END 
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In routine STGTST, and when ISTAGE=0, the routine may then be used to 
test for any of the increasing variables being greater than the values 
in the STEST array (i.e., increasing variable^ ^STESTi). In like manner, 
the decreasing variable^ > STESTDi test is made for i = 1, 2, 3, 4. 

Note that the testing stops at the first test to be satisfied. If a 
test is satisfied for increasing variables, the routine places the BCD 
name of the varible in STGVAR and the BCD word "INCR" in STGVAR-1. This 
process is similar for decreasing variables, except that "DECR" is 
placed in STGVAR-1. 

Similarly, when ISTAGE t 0, the routine may be used to test for al1 
of the increasing variables being equal to the values in the STEST 
array (i.e.,increasing variable i = STESTi). In like manner, the decreasing 
variable^ = STESTDi test is made for i = 1, 2, 3, 4. Note that in order 
to pass the test the conditions must be met for all i. 

c. Usage 

(1) Initialization - The statement, CALL STGTSI 

must be given. This statement must be given each time a new 
set of variable names is to be tested. 

(2) Testing Staging - To test if staging is indicated, 

CALL STGTST (INDSTG) 

The routine returns INDSTG = 0, No staging 

INDSTG = 1, staging indicated 
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FLOW DIAGRAM (STGTST) 
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4. INUPD, LNUPD, INPUZ, INTEG, UPDATE - INTERFACE ROUTINES FOR 
INTEGRATION ROUTINE 

a. Purpose - To serve as an interface between the integration routine 
proper (MIMIN) and any routine requesting a variable to be integrated. 


There are five logical functions which these routines perform. 

For a particular call to one of these routines, one of these functions 
will be enacted. The P array is the array of current derivative values 
of the variables that are being integrated. The Y array is the array 
of the current integrated variable values. 


b. Usage - for each of the interface routines. 

(1) CALL INUPD(N.L) 

The number of integrated variables is increased by N. The subscripts 
in the P and Y arrays for the values XDOT and X respectively are stored 
in the array L. 

(2) CALL LNUPD (M ) 

The number of integrated variables is decreased by M. 

(3) CALL INPUZ 

The P an -1 Y arrays are set to 0; the number of integrated variables 
is set to 0. 
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(4) CALL INTEG (K,XDOT) 

The value of XDOT is stored in P(K). 

(5) Call UPDAT (JX1, JX2, XJ1, XJ2, XJ3, XJ4, XJ5) 

JX1 = No. of variables (1 <_ JX1 £ 5); JX2 = subscript of first 
variable (XJ1); 2 nc * variable (XJ2) has subscript of JX2+1, etc. 

When the logical variable 3KIPUP (DIRCOM) is false, then the values 
of the integrated variables are picked up from the proper places in the 
Y array and put in XJ1, XJ2, etc. 

When SKIPUP is true, the values of XJ1, XJ2, .... XJ5 are stored 
in Y (JX2), Y (JX2+1), . . ., Y (JX2+JX1 - 1) 

INUPD will terminate the case if more variables are requested to be 
integrated than there is room for in the integration arrays; at the 
present, this upper limit is 50 variables. 

5. MIMIN - INTEGRATION ROUTINE 

a. Purpose - To perform the calculation necessary to integrate 
an arr$y of variables by the variable-step or fixed step Runge-Kutta 
Method; to determine an estimate of the relative error, and from this 
information, determine the new step size of integration in the 
variable-step mode of integration. 
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RETURN 


RETURN 













■■MailaaAtHdHMMHHtMMHIIIMli 
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b. Usage - The entry to this routine is as follows: 

CALL MIMIN (SN) 

where N = statement number to which nonstandard return is made 
from MIMIN. 

c. Method - The variable step Runge-Kutta Method is calculated as 
follows: 

STEP I 


Y mox n = | Y n | 

Y n = MX,Y n ) 

Y An = Y n ♦ h ' 2 Y n 
Y An s F (X ♦ h/2 , V An ) 

V B„ = Y n + “/2 Y An 

Y Bn = F < X + h/2 -W 
Y Cn = Y n ♦ h Y Bn 

Y Cn 2 F < X + h < Y Cn> 

Y Pn S Y n ♦ h ' 6 < Y n + 2 Y An ♦ 2 f ♦ Y- ) 


STEP II 


II 

c 

> Q 

Y n + h/4 Y n 

Y On- 

F(X + h/4,Y 0n ) 

II 

C 

Y n + h/4 Y D „ 

II 

c 

F (X + h/4 , Y Efl ) 

II 

c 

Y " + h/2 Y En 

II 

F(X ♦ h/2 , Yp n } 

Y ln = 

Y n + h/2 • 1/6 { Y, 


♦ Y F J 
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STEP III 


Y Gn 

s 

Y in + h/4 Y ln 

Y Gn 

= 

F(X + 3/4h, Y, ) 
on 

Y Hn 

= 

Y ln + h/4 Y Gn 

Y Hn 

= 

F(X + 3/4(i t Y Hn ) 

Y In 

= 

Y in ♦ Y Hn 

Y In 

= 

F ( X + h , Y In ) 

Yn 

= 

Y ln + h/2 - 1/6 <Y ln + 2 Y 6n ♦ 2 Y Hn ♦ Yj„) 


where n = 1, 2, 

3, . . ., NDEFEQ 



NDEFEQ = No. of 

differential 

equations 



Computation of 1 

Relative Error 

E R = ' 

iV |Y » - W 


Z n = Max 

(*-.„■ M 

) ; R = Mon 

[ Rt Hr 1/ Mo* ( 

z n. 1 )] 

d. Input for Integration 

Routine 



Symbol Used 

Math 

Symbol used 

Nominal 


by READ Routine 

Notation 

by Integration Value 

Remarks 

IVARBH 


IVARBH 

0 

Use variable-step 





=1, Use Fixed Step 

TIME 

t 

X 

0. 

Time to begin integ. 

DELTS 

At 

ox 

.1 

Time internal to int. 

AMINER 

At mi n 

DXMIN 

.001 

Minimum At 

AMAXER 

At max 

DXMAX 

10000. 

Maximum At 

RELER1 

^EL 

tL FRl 

RELER1 

.00007 

Rel. error tol. #1 

RELER2 

r el 

tL ER2 

RELER2 

.000005 

Rel. error tol. #2 


N 

N 


No. of diff. eqs. 


Y i 

Y(N) 


Array of dep. -var. 


P i 

P(N) 


Array of 

PRTMIN 


PRTMIN 


Print Minimum 



INDLG 


Landing gear indicator 



INDSTE 

1 




ST0PIT 





SW 
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FLOW DIAGRAM (MIMIN) 
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FLOW DIAGRAM (MIMIN) 
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6. LGDET - ROUTINE TO RESTRICT LG VARIABLES IN INTEGRATION ROUTINE 

a. Purpose 

When the statement CALL INTEG (K, XD) is executed, the derivative 
XD is stored in an array P at location K of that array (i.e., P(K)). The 
integration is not processed until a complete pass has been made through 
the program and all calls to INTEG have been made. 

As the integration is performed the Integrated variables are stored 
in an array Y at the correspondin' location K as its derivative (i.e.,Y (K)). 

When the statement CALL UPDAT (N,K,X) is executed, the integrated 
variable X Is transmitted from the Y array (location Y (K)) to the location 
X. The N designates the number of variables to be transmitted. 

Due to the requirements of the landing gear problem some variables 
that are integrated are restricted to certain conditions. Therefore, 
it was necessary to write the routine LGDET to restrict these variables 
in the integration routine. As stated above these variables are stored 
in the P and Y array in the integration routine. 


b. Linkage - CALL LGDET 
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FLOW DIAGRAM (LGDET) 
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FLOW DIAGRAM (L6DET. 




RETURN 
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7. ASRCH, TDATA - DIRECTORY SEARCH ROUTINES. 

a. Purpose - To provide a BCD word look-up from a subscript. 

b. Method - Given a subscript, the routine will search the directory 
for the BCD word corresponding to that subscript. If the subscripts do 
not compare the BCD name is set to blank and return to the Calling Program 
is made. 

c. Usage - Entry is made to the routine with the following statement 
a. CALL ASRCH (LOCI, SYM1) 

where 

LOCI = subscript being searched 

SYM1 = the variable name into which the routine is to store 
the corresponding BCD name. 

d. A call to the subroutine TDATA initializes table name data 
which is used by subroutine TABRE (in Overlay (1,0)). 
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8. DEF - HEADING AND PAGE EJECT ROUTINE 

a. Purpose - To provide page ejection and title printing. 

b. Method - Initially the current page number (NPAGE) is incremented 
by 1. The title is printed and returned to the calling program. 

c. Usage - Entry is made via the statement 
CALL DEF 

9. STFL, STFLD, STQVAR, ARRAY - Storage and Printing 
of Output Routines 

a. Purpose - To provide a method of printing output names of 

variables and thei* 1 values when necessary. Names of variables or values 
which are to be printed are not actually printed by the routines until 
at least eight have been accumulated by a series of calling sequences. 

b. Method - Each time a call is made, names of variables or values 

are saved until eight are stored. At this time they are printed. This 
process is repeated until all names or values have been handled. If 
less than 8 remain, they are saved for further call statements or until 
forced to be printed. 

c. Usage 

(1) The printing of Hollerith Code 
CALL STFL (JOPT, N, ARG1) 

JOPT = 0: Force any possible remaining print 

JOPT = 1: See (2) below for printing of Values of Variables. 

JOPT = 2: Prints N words of Hollerith information from ARG1 (1), 

ARG1 (2), . . ., AP.G1 (N) where there is at most 6 characters per word. 
JOPT = 3: Prints 1 word of Hollerith information. 

N should equal 1. 
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The above call adds the N Hollerith code words to the list of code 
words to be output on the next line of print. When 8 code words have been 
accumulated, the line is printed; any excess code words are added to 
the list for the next line of print. Codes are printed with "7X,A6,7(9X,A6)" 
format. 


(2) The printing of Values of Variables 

CALL STOVAR (N, A1, A?, . . ., A8) 1 < N < 8 
where N is an integer identifying the number of arguments following it. 
If N <_ 8, then the remaining arguments must be dummy arguments. For 
example, CALL STOVAR (6, A1, A2, A3, A4, A5, A6, DU, DU). 

CALL STFL (JOPT, N, ARG1). This call prints the value 
of one variable when JOPT = 1, N = 1, and ARG1 is the name of the 
variable. 

Lines accounting is taken care of within this routine. 
Values are printed with a "1PE15.7" format. 

(3) Forcing final print 

To force any possible remaining print 
CALL STFL (0, 1, DU) 
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FLOW DIAGRAM (DEF.) 



NPAGE « NPAGE+I 
LONG « t 
EJECT PAGE 

write: “six-degrees 

OF FREEDOM FLIGHT PATH 
STUDY GENERALIZED 
COMPUTER PPOGRAM" 

WRITE " INDSOF Z" 
CASE (NCASE) ■STAGS 
INSTAGE) . PAGE (NPAGE) 
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' STFL > 
, ENTRY j 


J#PT 

I *O 


GETARG, • ARG, 
WENT * 5 
K • 2 

CALL ARRAY<N,0) 


RETURN 


CLOW DIAGRAM (STFL) 


/STFLD\ 

[entry/ 


INTIALIZES BY 
DATA STATEMENT: 
IZ * O 

CLEAN *.TRUE. 
INTEG = FALSE. 


ARG 2 

ARG 1 j 

... N 





{ *91 > 


NENT * I 


NENT » 2 

/ N 

1 

' 

K * 1 

- ^return) 

CLEAN 

* TRUE. 1 

CALL ARRAY(N.O) 


CALL ARRAY ( N, 1 > | 



(RETURN 


RETURN Hs-- 
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FLOW DIAGRAM (ST$VAR) 
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FLOW DIAGRAM (ARRAY) 



b2 















FLOW DIAGRAM (ARRAY) 



ALIST (12 ) * GETARG (J) 
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10. LINES - LINES ACCOUNTING ROUTINE 

a. Purpose - To keep an accounting of the number of lines printed 
per page, and to provide for page control. 

b. Method - If the number of lines to be printed (LC0UNT) is such 
that it will not fit on the current pane, the page is ejected (via DEF) 
and printing will begin on the new page. Initially, the location L0NG 
should be set to zero, indicating that currently no lines have been 
printed on the present page. 

c. Usage - Entry is made via the statement, 

CALL LINES (LC0UNT) 

where 

LC0UNT = A fixed point variable or constant indicating the 
number of lines to be printed. 

11. ASIN - ARC SINE FUNCTION 

a. Purpose - To compute the arc sine of a normalized floating 
point argument (X). 

b. Method - , , 

Sin" 1 (X) = J - Cos -1 (X) 

c. Usage - The Arc Sine is computed using the statement, Y = ASIN (X) 

where jX| < 1., and Y = Sin"^ (X). 

12. ACOS - ARC COSINE FUNCTION 

a. Purpose - To compute the arc cosine of a normalized floating point 
argument X. 

b. Method - For |X| < 7.4505806 x 10 the arc cosine Is set equal 
to tt/ 2. For X = 1., arc cosine is set equal to zero, and for X = -1 arc 
cosine is set equal to tt . When the argument X / +1., the routine gives 
the arc cosine in radians from C to tt. 
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FLOW DIAGRAM (LINES) 
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FLOW DIAGRAM (ASIN.) 
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c. Usage - The arc cosine is computed using the statement 
Y = AC0S (X) 

where |X| £ 1. and Y ~ Cos" 1 (X). 

13. ATAN2 - ARC TANGENT ROUTINE 

a. Purpose - To compute the arctangent of the quotient of two 
normalized floating point quantities Y/X, with proper quadrant control. 

b. Method - The routine computes the quotient Y/X. The arctangent 

is computed with quadrant according to the sign of Y and X. If X = 0 

-1 Y * 

and Y f 0, the routine computes 1 - Tan (Y/X) = j y j* . 

If X = 0 and Y = C, it computes 
Z = Tan" 1 (Y/X) = 0. 

c. Usage - The arctangent of Y/X is obtained via 
Z=ATAN2(Y,X) 

14. ERROR, EXERR - GENERAL TABLE ERROR ROUTINE 

a. Purpose - To provide a method of indicating the table which 
may possibly contain an error. Also, to provide the Stop Number that 
identifies an error condition in an equation. 

b. Method - By the statement CALL ERROR (LOCT) in which LOCT is 

the subscript of the curve in error, the routine will search the subscript 
table and find the corresponding BCD word. This word will then be 
printed as: 

"TABLE ERROR AAAAAA" 
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Where AAAAAA is the BCD name of the table. If the name cannot be found 
in the directory 

"TABLE ERROR 

LOCATION OF TABLE NOT LISTED IN DIRECTORY" is printed and a 
return to the calling program is made. In either case INDSTE is set 
to zero. By the statement CALL EXERR (NUM) in which NUM is the stop 
number, the routine will write "STOP NUMBER III" 

where III = the stop number. If NUM = 0, an exit is made from the routine 
with no printing. In either case INDSTE is set to zero. 

c. Usage - Entries are made to the routine with the following 
statement: 

(1) CALL ERROR (LOCT) 

where LOCT is the table subscript. 

(2) CALL EXERR (NUM) 

where NUM is the stop number. 

15. NDTLU - N-DIMENSIONAL TABLE LOOK-UP ROUTINE 

a. Purpose - To provide a method of linearly interpolating in a 
table of n independent variables. 

b. Method - Given the arguments X(l), X(2), .... X(N-l), the 

routine computes Y = F(X(1), X(2).X(N-l)) by linear interpolation 

from a table. 
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c. Usage - Entry is made via the statement: 

CALL NDTLU (ND, NA, X, Z, XA, ZR, IE) 

where 

ND = Dimension of lock-up (when V = F(X), ND = 2) 

NA - An array of length ND-1. Numbers of values of each table of X. 

The tables are listed by size. 

X = Tables of each X in order. 


Z = Functi 

on Vali 

ues. 

If 

array 

must be ir. 

the 

Assume NX=4, NY 

=3, NZ ; 

=2. 


W(l) = 

F(XI, 

Yl, 

Zl) 

W(2) = 

F(X2, 

Yl, 

zi) 

W(3) = 

F(X3, 

Yl, 

Zl) 

W(4) = 

F(X4, 

Yl, 

zi) 

W(5) = 

F(XI, 

Y2, 

Zl) 

W(6) = 

F(X2, 

Y2, 

zi) 

W(7) = 

F(X3, 

Y2, 

Zl) 

W(8) = 

F(X4, 

Y2, 

Zl) 

W(9) = 

F(XI, 

Y3, 

Zl) 

W(10)= 

F(X2, 

Y3, 

Zl) 

W(ll)= 

F(X3, 

Y3, 

Zl) 

W(12)= 

F(X4, 

Y3, 

Zl) 


ZR = Results 
IE = Error Code 
= 0 No error 
-1 X array toe small 
1 X array too large 


F(X, Y, Z) the dependent variable 
owing order. 

W(13)= F(XI, Yl, Z2) 

W(14)= F(X2, Yl, 12) 

W(15)= F(X3, Yl, 12) 

W(16)= F(X4, Yl, 12) 

W(17)= F(X1, Y2, 12) 

W(18)= F(X2, Y2, Z2) 

W(19)= F(X3, Y2, Z2) 

W(20)= F(X4, Y2, 12) 

W(21 )= F(XI, Y3, Z2) 

W(22)= F(X2, Y3, 12) 

W(23)= F(X3, Y3, U) 

W(24)= F(X4, Y3, 7.2) 


2 array not ir. ascending order 
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Let n be number of indpendent variables, then the table is called an "(n+1) 
dimensional table." 

Z = F(X 1 , . . . , X n ) 

To use a table of dimension 1 3 ana ± 5, a call to HIH0 should be 
made with the list of arguments in the calling sequence in the same order 
as the independent Variables are numbered. 

16. ATMS - ATMOSPHERE CALCULATION ROUTINE (1969) 

a. Purpose - To compute the atmosphere characteristics: Density, 
speed of sound, pressure, temperature, and kinematic viscosity. All 
are a function of altitude. 

b. Method - All atmosphere characteristics are computed using the 
1969 ARDC model atmosphere. Values of the atmosphere characteristics 
are computed for positive altitudes. If an altitude is negative, the 
sea level value will be obtained. 

c. Usage - Linkage is affected by 
CALL ATMS (HGC7F) 

where HGC7F = altitude in feet 

17. INVR3 - INVERSE OF A NONSINGULAR 3X3 MATRIX 

a. Purpose - To compute the inverse of a nonsingular 3x3 matrix. 

b. Method - Let a = [a..] i, j = 1,2,3 

• J 

then A"'* - ( ii [A..] is computed where A., is the cofactor of a... 
The matrix A must be stored in the usual FORTRAN sense (i.e., 
columnwise). 
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FLOW DIAGRAM (NDTLU) 



RETUR 
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MB (1) * 

0 . 

MB( 2 ) - 

11000 . 

MB(3) ■ 

25000 

HB(4) = 

47000. 

MB ( 5) * 

53000. 

HB ( 6 ) = 

79000. 

HB ( 7) = 

90000. 

MB( 8 ) * 

105000. 

HB (9 ) * 

160000. 

HB (10) = 

170000 

HB (It) ' 

200000 . 

RH^B(l) 

= 23.7692E - 4 

RH$B(2) 

= 7.0620 E -4 

RH0B(3) 

= 7.7650E -5 

RHWB(4) 

= 2.8804E-6 

RH 9 B(5) 

= I3.9468E - 7 

RH$B( 6 ) 

= 4.II89E -8 

RH$B(7) 

= 4.26I0E-9 

RH$B( 8 ) 

= 2.2320E -10 

RH$B(9) 

= 1 8450E-I2 

RH()B (10) 

= 1 .3380E-I2 

RH(^B (II) 

= 6.II30E-I3 

PB(1) = 

2116.21695 

PB(2) = 

472 .73 

PB(3) = 

51.979 

PB(4) = 

2.5155 

PB(5) * 

1.2181 

PB ( 6 ) = 

2.I08E -2 

P8(7) = 

2I.809E -4 

PB( 8 ) * 

I5.562E -5 

PB(9) = 

75.578E-7 

PB(IO) = 

58.954E -7 

PB (II) - 

29.759E-7 

TB(1) = 

518.688 

TB(2) = 

389.988 

TB( 3) * 

389.988 

TB{ 4) = 

LOS.788 

T 8 (5) = 

508.188 

TB( 6 ) = 

298.188 

TB( 7) * 

298.188 

TB ( 8) * 

406,'I88 

TB( 9) = 2386.188 

TB(IO) =2566.188 

TB (II) =2836.188 


FLOW DIAGRAM (ATMS) 


AK1 (1) = 

- 225569E-4 

AKi (2) = 

.0 

AKi (3) = 

.138466 E-4 

AKi (4) = 

.0 

AKI (5) = 

-.159202 E-4 

AKI (6) = 

.0 

AKi (7) = 

.24I458E -4 

AKI (8) = 
AKI ( 9) = 

. 886289E-4 
.T54341E-5 

AKI (10) = 

.350715 E-5 

AKI (II) = 

. 222I29E -5 

AK 2 (1) = 

- 5.25612 

AK 2 ( 2) = 

0. 

AK 2 (3) = 

11.3883 

ii n 

ift 

NM 
K X 
« < 

0. 

- 7.59218 

AK 2(6) = 

0. 

AK 2 (7 ) = 

8.5412 

AK 2 (8) = 

1.70824 

AK2(9) = 
AK 2 (10) = 

3.41648 

6.83296 

AK 2 (11) = 

9.76137 

AK 3 ( 1) = 

.0 

AK 3 ( 2) = 

I57689E-3 

AK3 (3) * 

.0 

AK 3(4) = 

.I20869E- 3 

AK3(51= 

.0 

AK 3(6) = 

.206234E-3 

AK3(7)= 

.0 

AK3 (8) = 

.0 

AK3 (9) = 

.0 

AK3 (Id = 

.0 

AK3 (II) = 

.0 

All) 

.759511 

A (2) = 

.935787 

8(1) = 

.174164 

B (2) = 

.273966 

C(i) = 

220. 

C (2) = 

180. 

D(i) = 

25. 

D (2) = 

140. 


CO 
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INDATM * 0 


.MGC7F > 2.5E6. 


FLOW DIAGRAM (ATMS) 


HGP > 90000, 


VS77F * 0. 
TA77R = 0. 
PA77P = 0. 
ANUA7F = 0. 
RH$AS = 0. 


RETURN 


HGT90 = . FALSE. 


HGP < HB (11) 


LAY = 11 


HGC7F < 0. 


VS77F = 1116.43372 
TA77R = TB(1) 

PA77P = PB(1) 

ANUA7F = I.3723288E-4 
RH^AS = RH$B (1) 


1=1 + 1 


,^-(HGP>HB ( I) 


LAY = I -1 


TMP = .3048* HGC7F 
HGP = TMP/I. + 
TMP/6356766.) 

HCT90 = TP'.'E 
m = a 


HGP > 180000. 


i a m 


TMP = HGP - HB (LAY) 
TMP2 = I.+ AK1 (LAY * TMP 
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FLOW DIAGRAM (ATMS) 
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c. Usage - Linkage is obtained via the statement: 

CALL INVR3 (A,B, INDER) 

where 

A = the array name of the matrix to be inverted. 

B = the array name where the resulting matrix is to be stored. 

INDER = an error indicator set by the routine. 

(a) INDER = 1, results are good 

(b) INDER = 2, A was 0. 

18. MULT31 - A MATRIX MULTIPLICATION ROUTINE 

a. Purpose - To postmultiply a 3 x 3 matrix by a 3 x 1 matrix. 

b. Method - The result of [A] [B] = [C] is computed using single 
precision floating point arithmetic. All elements must be stored in 
the normal FORTRAN sense (i.e., columnwise). 

c. Usage - The matrix multiplication is obtained by the statement: 

CALL MULT31(A,B,C) 

where 

A = array name of the 3x3 matrix [A] 

B = array name of the 3 x 1 matrix [B] 

C = array name of the resulting 3 x l matrix [C] 


71 
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FLOW DIAGRAM (IjjjVRS) 
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FLOW DIAGRAM (MULT 31) 
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19. TRNP0S - A 3 x 3 MATRIX TRANSPOSE ROUTINE 

a. Purpose - To transpose a 3 x 3 matrix [A] to obtain he 3x3 
matrix [A]'. 

b. Method - The resulting transposed matrix is stored in a 
separate array. All elements of [A] must be stored in the normal 
FORTRAN sense (i.e., columnwise). 

c. Usage - The transpose of a 3 x 3 matrix [A] is obtained by: 

CALL TRNP0S (A,B) 

where 

A = The array name of the 3x3 matrix [A]. 

B = The array name of the 3x3 matrix [A]. 

20. HIH0 - N-DIMENSIONAL TABLE CALL ROUTINE 

a. Purpose - To set up the NA array and Z location of tables with 
dimension from 3 to 5 as required by the calling sequence to NDTLU, 
which is 

CALL NDTLU(ND,NA,X,A,XA,ZR,IE), to make the call to NDTLU, and return 
the function value or data on a table read error. 
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FLOW OtACRAM (TRNP4>S) 
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b. Usaqe - Linkage to the subroutine is made via the statement 
CALL HIH0 (N, L0CT, NX1, NX2, MX3, NX4, X1ARG, X2ARG, X3ARG, X4ARG, A) 
where 

N = dimension of table look-up, when A = F(x), N=2. 

L0CT = location of the first value in the table. 

NX1 to NX4 - location cf number of points in the XI to X4 
array of independent variable values 
X1ARG to X4ARG - name of XI to X4 argument or a dummy location if N < 5 
A = location of the dependent variable. 

21. TLU - TWO-DIMENSIONAL TABLE LOOK-UP ROUTINE 

a. Purpose - Given an argument X, to compute Y = F(X) from a table 
of X and Y value: <v linear interpolation. 

b. Method - The table of X values is searched until for some 
i, X ( < X < X u , 

Linear interpolation is then performed. If, for some i, X=X i then Y 
is set to Y^. 

c. Usage - Ent^y is made via the statement, 

CALL TLU (X, L0CT, Y) 

where 

X = = variable name of the argument 

L0CT = location of first subscript of desired table of data 
Y = variable name of the interpolated value. 

All tables are stored in a table called C. Therefore, i.he first value 
of a particular table is located at C(L0CT). As an example, assume 
foat the name of the desired table is FTAB01. Then the curve must 
be stored as follows: 

C(L0CT) = FTAB01 (1) - N=iio. of points in curve (integer) 








AFFOL-TR-71-155 
PART IV 
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FLOW DIAGRAM (HiHifr) 


( HIH(J 1 
iENTRY , 


XA (1) 

= 

XI AR6 

XA( 2) 

s 

X2 ARC 

XA(3) 

s 

X3 ARC 

XA<4) 

= 

X4 ARG 

NA (i ) 


NX 1 

NA(2) 

= 

NX 2 

NA'35 

= 

NX 3 

NA(4) 

= 

NX4 

NAT 

* 

NX 1 


IERR$R = O 

Tno 


CALL^IL 
CALL ERROR (LfrCT) 
CALL LINES (2) 
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C(L0CT+1) = FTAB01 (2) = X ] 
C(L0CT+2) = FTABO'I (3) = Y ] 
C(L0CT+3) = FTABOl (4) = X 2 
C(L0CT+4) = FTABOl (5) = Y ? 


C(L0CT+2N) = FTABOl (2N) = X N 
C(L0CT+2N+1)= FTABOl (2N+1) « Y^ 

22. TFFS1 - ENGINE THRUST AND THROTTLE SETTING 

a. Purpose - To provide a method of introducing the engine 
thrust characteristics into the computation, with an option to find the 
throttle setting that corresponds with a certain thrust and Mach number. 

b. Usage - Linkage to TFFS is accomplished via the following 
statements: 

(1) CALL TFFS ! 

Pre-data initialization. 

(2) CALL TFFS3 

Thrust computation section. 

(3) CALL TFFS4 
Initial print 

(4) CALL TFFS5 

Code printing to identify the coming time history. 

(5) CALL TFFS6 

Time History Print 
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(6) CALL TFFS7 

Update integration (none for this subprogram) 

23. VPCS - Vehicle Physical Characteristics 

a. Purpose - To introduce various physical characteristics into 
the general solution of the problem. Here, mass, moments, and products 
of inertia of the vehicle and rotating machinery, reference lengths aid 
areas for aerodynamic coefficients, jet damping characteristics lengtns, 
and center of gravity information are included. 

b. Usage - Linkage to VPCS is provided ly the following statements: 

(1) CALL VPCS1 

Pre-data initialization. Necessary romina values are set. 

(2) CALL VPCS2 

Post-data initialization. 

(3) CALL VPCS3 

VPCS computations are performed if INDVPC is non-zero. 


BO 


*4 
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I 



FLOW DIAGRAM (TFFS1) 



SET TO zero: T, , 

Ty , Tj , Ly , My t Ny 
TXB7P = 0. 
TYB7P = 0. 
TZB7P = 0. 
ALT77F = 0. 
AMT77F= 0. 
ANT77F= O. 


RETURN 




SET TO ZERO: T, , 

Ty i Tj | Ly | My t Ny 

TXB7P= 0. 
TYB7P= 0. 
TZB7P * 0. 
ALT77F * 0. 
AMT77F = 0. 
ANT77F * 0. 


™_(„,o777o") 

| N0 

I = 0 


1*1 + 1 

OBTAIN T(I)-/(N(I),M n ): 

CALL HIH<t> 

MT(I) ~ T(I) » ZN ( I) 

NT (I) = -T(I)* YN (I) 
TXB7P = TXB7P + T (I) 
AMT77F * AMT77F + MT{I) 
ANT77F - ANT77F +NT( I) 



} 


i 

{ 

\ 


bl 


Mir-- ■ i 
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FLOW DIAGRAM (TFFSI I 
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AFFDL-TR-71-i5B 
PART IV 


FLOW DIAGRAM (VPCSI) 



SET TO ZERO ! I xx , 
’yy ■ 1 1* > *xy ! yz > 

i •^mi^ir'^y’^Z 
BSIIXX = 0 . 
BSIIYY = 0 . 
BSIIZZ =0 

bsitxz = o 

BSIIXY S 0 . 
BSIIYZ* 0 . 
ALLJDF= 0 . 
AIMJOF-O 
ALNJDFsO. 
ALYJDF =0 
A 1 _ZJDF= 0 . 


Q INOXZS-O 


YES 


NO 



A1XY 

AIYZ 

BS-0 
BS = 0. 
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FLOW DIAGRAM (VPCSI ) 
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(4) CALL VPCS4 
Initial print 

(5) CALL VPCS5 

Code printing to identify the coming time history is performed 
if INDVPC f 0. 

(6) CALL VFV6 

Time history print. 

24. SACS - AERODYNAMIC FORCES AND MOMENTS 

a. Purpose - To provide a complete accounting of the various 
contributions to the aerodynamic forces and moments, regardless of the 
flight conditions or the vehicle considered. In SACS coefficients 
are computed in the proper coordinate system for use in other parts of 
the TOLA program. 

b. Usage - Linkage to SACS is accomplished via the following 
statements: 

(1) CALL SACS1 

Pre-data initialization. 

(2) CALL SACS3 
Aerodynamic computation. 

(3) CALL SACS4 

Initial print. (None for this subprogram). 

(4) CALL SACS5 

Code print to identify the trajectory point for SACS. 

(5) CALL SACS6 

Time history print for SACS. 

(6) CALL SACS7 

Update integration (None for this subprogram). 
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(7) CALL SACS8 (C LR , a , C DR ) 

Determine the a d that corresponds with and Cp. 

(8) CALL SACS9 (o , C L> C D ) 

Determine the C|_ and C Q that corresponds with a . 

(9) CALL SACS10 (a ) 

Determine 6 as a function of a . and all current variables, 
qri d 

(10) CALL SACS11 (<$ rN ) 

Determine 6 „ as a function of all current variables. 
r N 

25. AER01 - AERODYNAMIC DATA LOOKUP FUNCTION 

a. Purpose - To look up aerodynamic data from the table array in 
C0MMON TBDIR/C(300). It looks up the first two values of a particular 
table depending on the argument of the function. 

b. Usage - Linkage to AER01 is accomplished in the following ways: 

Y = AER01 (L0CT, AER02) 

where L0CT = subscript for the table C(300) 
which will be the 1st location of a particular table that one is interested 
in. 

After the function is executed, 

AER01 = 1st value of the table that begins at C(L0CT) 

AER02 = 2nd value of the table that is located at C(L0CT + 1) 

26 . AQUAD - A QUADRATIC FUNCTION 

a. Purpose - To solve the equation a{x|x + b x + C = 0 

b. Usage - Linkage to AQUAD is accomplished in the following way: 

Y = AQUAD (A, B, C, XLLIM, XULIM) 
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SACSI 

ENTRY 


SET 7® ZERO. a,n F , 

y ■ £ t w i n 

AA77P= o. 

AN77P - 0. 

YA77P = 0 
ALA77F = 0. 
AMA77F a0 
ANAZ7F = 0. 


- 

... 

URN^. 

YES 

wu r N 

1 NDAER= 0 

] 



ENTRY 

SACS4 


ENTRY 

SACSS 


FLOW DIAGRAM ( SACS I) 


ENTRY 
SACS 7 


WRITE heading;"SAC l' 


INDAa0=0 J 

NO 

WRITE C®OE! “CAVAH” 


WRITE CDDES; “CA 
CN CY" 

WR.TE COOES.' "CL 
CM CNN" 



RETURN! 


® WRITE VALUES FOR! CA, 
CN, CY 

• WRITE VALUES FDR; CRM, 
CM,CYM 


WRITE VALUE FOR;CAVAH 

NO 

INDAao=Q 

NO 

(~ IN0AER = 0 




8i 
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FLOW DIAGRAM (SACS: ) 


(SACS9) 

* 3 


0 


YES 


/-\ NO 

( INDAOI = 0 )—*- C A = (AER®I(L©C(I),AER®2)-AER®2)»TMRI 
-|- ' ___+AER®2 

c 


INDA02- 0 


YES 


C A a * = (AERCI (LCC<2), AERC2) ~ AER®2)*Tw»P<l) 

+AER®2 


INDAI 

YES 

l 

5=0 ^ 

1- 

NO 

Cm s =( AERCI (L®C(I5),AER®2)-AERC3*TMP(I) 
9 +AER®2 

b 

__ 1 

— - 

[-- N0 1 

1- -1 

i 


YES 


/-*- 

( INDAI6 = 0 V-•> C NJ s 2 = <AER®I(LCC(I6), AER®2)-AER®2)*TMP\I) 

V--- S oq +AERC2 


YES 


/ --NO 

( INDAIO-O V-C N = ' AERCI (L©C( 10), AERC2)-AER«2)*TMP(I) 
V - ' +AERC2 


I NDAI I =0 ^ 


NO 


YES 


c 


f 


C N = (AEROI (L®C( I I), AER®2) - AER®2)-*-TMP(l) 

+AERC2 


INDAi 2-0 


NO 


YES 


C N ~.AERCI (L®C( 12), AER®2) -AER®2 )*TMP(I) 
a +AER®2 


c 


NO 


INDA20 =0 


YES 


C N * =(AER®i (L®C(20), AERC'2) - AERC2)* TMPO j-. 

+AERC2 1 


-L_- NO 

^ INDA22=0 J —►- 


f"- 

/— 1 —\ r 

( INDA80 =0 )—fc-J 


C. N - ! AER® I (L®Cf 22). AER®2)- AFR®2)»TMF<I) 
** +AERC2 


c A n = (AERCI (L®C (70), AER®2) -AERC2)»TMP(I) 

+AERC2 


0 


Y -^G) (SA0S8I 
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FLOW DIAGRAM (SACS!) 
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FLOW DIAGRAM (SACSI) 


9 


r 


y~® 




1602 


(SACS10) 


r - 

INDA2 

V- 

YES ] 


NO 

C y £ = (AER®I(L®C(27),AER®2)-AER«2)*TMP(I)+AER®2 

1 

□ 


-. NO 1 

1 - -1 

f 


'NDA28=0. 

"vcsl 


3' 


Cy£*= (AER«I (L®C(28) ; 4ER<r2)-AER®2)*TMP(l)+AER®2 


1 


S~ -*- n no 

( INDA68= 0 V-*. 


YES 


C„p ~! A£R®I (LOCI68), AERC2) - AER®2)*TMP;i)+AER®2 


□ 


Q INDA< 

-- NO 

>9 r 0 V— *»- 

C, pH - (AER®I (L®C(69), AEP02) - AER®2)*TMP< 1 )+AER®2 

— 

YES 



_J 

--3 

-^ NO 

1— 

1 


V. 


YES 


r 


E 


Nl 


V 


INDA7I = 0 


YES 


C n g r 2 = (AER®I (L®C( 71), AER®2) - AER®2)*TMP(I) +AER®2 


( I GO 3 


(SACSII) 


= 1 


- 1 -.NO ! - 

^ IN0A03= 0 J — C/^p = (AER®I(LOCO), AER®2)-AER®2)*TMP(4)+AER©2 

YES I '—-- 




1 


NO 


IN0A04 -0 


Ca^ 3 ! = 'AEROI (l®C(4).AER®£)-AER®2!*-TVP(l)+AER®2 
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FLOW DIAGRAM (SACSI ) 
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FLOW DIAGRAM (SACS I) 
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FLOW DIAGRAM (SACS I) 


INDA34 = 0 
YES I 

INDA35 z O 

YES ^J 

IN0A36=O 
YES I 

INDA37 = 0 

YES ^— 

INDA38=0 

YES ^*_ 
INDA39=0 
YES ^ 

INDA40=0 

INDA4I =0 
YES m 

INDA42 - 0 

yesT"~ 


C y 0 = (AEROi (L®C(34), AER®2) -AER®2) *TMP(I)+AER©2 

■ Cy^ = (AER®l(L®C(35),AER®2)-AER®2)*TMP(l)+AER®2 

C yr = (AEROI (L0C(36), AER®2) -AER®2)*TMP(I) +AER®2 
_ 

Cy rx =(AER®i(L®C(37), AER®2)-AER®2)vTMP(l)+AER®2 

C/ 0 - (AERCI (L®C(3B), AER®2) - AER®2)*TMP( I) +AER®2 
C la = (AER®I(L®C(39) < AER®2)-AER«2)»TMP(I)+AF.R®2 j - 
c; a * = (AEROI (L®C(40), AER®2)-A£R®2)*TMP( I) +AER02 - 

C/0 - < AER®I (L®C(4I), AER®2)-AER®2)*TMP(I) +AER®2 
C/0* = (AER®K L J® C (42),AER©2)- A ER®2)*TMP( | )+ A ER®2 - 


/- 1 -N 

{ INDA43 - 0 

^ -Ip ' 

NO 

-, 

C/Sp * (AEROI (L®C(43), AER®2)-AER®2 )*TMP(I)+AER®2 1 

n 

YES 1 _ 


-1 

r 



^ INDA44 - 0 J 

NO 

C /Sp* -(AER®I (L®C(44), AER®2) -AER®2)*TMP( U+AER02 
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AfFDL-TR-71-155 
PART IV 


FLOW DIAGRAM (SACSl) 


, ->--sNO 

7 INDA45 - 0 /-*- 


YES 

C INDA46 


0/g a = (AER®! (L®C (45), AER®2) - AER®2 )*TMP(I )+A£R©2 


YES 


c 


- N NO 

*° f-*- c /aSp" 


= (AER®I (L®C(46),AER©2) -AER®2)*TMP(I) +AER®2 


IN0A47 = 0 




NO 


YES 


C(L ® C(4T) > AER ® 2J ~ AER ® 2) * TMP<I)+AER ® 2 


<*- 

INDA4 

1 - N NO 

ie = o *- 

C/ = (AERO! {L©C(48) f AER©2)-AER©2)*TMP(I)+AER©2 

— 

YES 






-. NO 1 


i 


^ INDA49- 0 y 


YES 


c 


C/ f = (AEROI (LfflC(49), AER®2) —AER®2 )-*TMP(l) +AER02 


NO f 


INDA50* 0 


YES 


K 


C i = (AER0I (L0C(5O), AER02) -AER©2)*TMP( I) +AER02 
r x 


h 


INDA54 


YES 


c 


- N NO 

i °y~ 


C m/ g = (AF.R®I { L®C(54), AER®2)-AER®2)*TMP(I) +-AER®2 


NO 


INOA55 = 0 


YES 


C 


C m; g2=(AER®l(L®C(55), AER®2)-AER®2)*TMP(l)+AER®2j- —. 


INDA64 


^rv: 


NO 




YES 


C mgx -(AER®l(l.®0(64),AFR®2)-AFR®2)*TMPil)+AER®2 


^ INDA65=~Q Y 

t ' 


C no = (AER©l(UDC(Si>), ALR®2)-AER®2)«TMP(I)+AER®2 


^~INDA66^ 0 ^—»- | C nc 


: (AER®I (LOClG6), AER02) - AER02P*TMP( I) +AERf 
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FLOW DIAGRAM (SACSI) 
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FLOW DIAGRAM {SACSI) 



Gy = C yo + Cy a lal + C ya 2 a 2 + Cy£/3 ♦ C y/ g2/3l/3l 

♦ C y g r Sr + C y g f 2 |Sr) Sr + C ya g r |a| Sr + C ia fi\a\ (3 

+ C vv38r 101 Sr +{C iP + c y / 9 x ^ x C . G .) 0cl 2 / 2V A 

+ (C yr + C yry AXc Q ) r d 2 / 2 V A 

Crm =G/ 0 + C* a |aj + C* a 2 a 2 + j/3| + C/g p Sp 

+c i8p 2 N Sp + C* agp |a| Sp + C W3 l a 0| 

+ c /)3Spi0! 8p + C/ p Od 2 /2V A + (C/ r +C/ r)t AXc. G .)r0 2 /2V A 

Cm s Cm 0 "*■ C(n a Ct + Gm<j 2 a l a ! + ^m^Ql 01 + @ 

+ C m Sq S< » H C mSq 2 8< » l 8p l + C maSq i a l 8< * 

^" C rrva)3 “ 101 f C m£Sq l^l 8 * + ( C md + C md, AX Ce ) ad,/2V A 
+ <C m q + C mq>i AX C q )qd,/2V A ~C N AX C G /d. 


c ym = c r<o < -r.a l a l + c Ma 2 a * + C n/ g/3 + q^g £ + C n/ g2 0|0| 
+c n8r ^ + C n g r 2 !8r{ Sr ^C na g r !«ll Sr + C na/ Q|a|0 
+C n/38r I0l 8r ^ « c n/S +c n^x^ x C.G.)/3d2/2V A 
+(C nr + C nrx AX C G ) rd 2 / 2V A —Cy AX C G ./ d 2 
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FLOW DIAGRAM (SACSI) 


IAP < 3 ] 

" --- ' 

NO 

|Hg-H R |< 10“*) 

' [no 


ENTRY 
SAC59 , 


IAP< 3 

| NO ' 

|Hg-H R |<IO" # 

I NO 


1 G(D 1 =2 

-—_-_ 


- 1 _ 




1 G<D 1 = 3 
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FLOW DIAGRAM (SACSI) 




fil 
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FLOW DIAGRAM (SACS!) 



TMP2 » C n g r , q*Sd 2 
TMP3 • C ngr q*Sd 2 
TMP4 • [c„ p -r y ^ AX CG /d 2 

+ l/SdKCoyg* -C yp * AX CG /d 2 l] q*Sd 2 


SOLVE FOR Sr N WHEKE Sr L < Sr N S8r u: 
TMP2|8r N |8f N + TMP38r N +TMP4*0 
(THIS IS DONE BY USING THE FUNCTION AQUAD) 


-{RETURN 

v y 
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FLOW DIAGRAM ( AERO I) 
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where A = the coefficient of JX)X 
B = the coefficient of X 
C = the constant term 
XLLIM = the lower limit of X 
XULIM = the upper 1imit of X 

After the function is executived, the value of X will be stored in AQUAD 

27. OPTI - Six-Degree -Of-Freedom Trajectory Program Over a Flat Planet 
a. Purpose - This program permits the computation of the six rigid- 
body degrees-of-freedom of a flight vehicle where the motion is assumed 
to occur over a limited portion of the planet. Under this assumption, 
the motion is taken to occur within a rectangular coordinate system. 

This program of computation has been provided to aid in such analyses as: 

(1) Boost-phase dynamics. 

(2) Pitch-roll-yaw coupling motion investigations. 

(3) Autopilot response to parametric disturbance. 

(4) Landing, approach, and flare-out maneuvers, etc. 

(5) Landing gear dynamics. 

The effects upon the motion of the body of the following items may also 
be accounted for: 

(1) Three components of wind velocity and acceleration measured 
in a local-geocentric Cartesian coordinate system. 

(2) An arbitrary planetary atmosphere. 

(3) The gyroscopic torques imposed by rotating machinery aboard 
the vehicle. 
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HX,) = X llim (A|X llim | + B) +C 
F(X t ) * X uL j M (A|X UL1M j + B)+ C 
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(4) The action of an arbitrary* three-plane autopilot upon the 
control deflections of the flight vehicle in response to the measured 
motion. The motion may be obtained from a stabilized platform. In this 
program the motion is more conveniently printed out in the Cartesian 
Coordinates, 
b. Linkage 

(1) CALL 0PT1 

Pre-data initialization. 

(2) CALL 0PT2 

Initialization after data is read in and computation for 
initial time point. 

(3) CALL 0PT4 

Calculation of equation of motion, and all other variables 
that are desired. 

(4) CALL 0PT6 

Time history printout. 

(5) CALL 0PT7 

Update variables that are integrated by the integration routine. 

28. LGEARi - LANDING GEAR CALCULATIONS, PART I 

a. Purpose - The LGEARI subroutine computes the effects of ground 
reaction and landing gear dynamics on the motion of a landing vehicle. 

A maximum of five independent landing gears may be used. 

b. Usage - Linkage to LGEARI is provided by the following statements: 
(1) CALL LGEARI 

Pre-data initialization. 

* Currently, only a Pitch-Yaw-Roll sequence is programmed, due to 
storage limitations. 
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FLOW DIAGRAM ( $PTI ) 


4. 


! 


I 



3 
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FLOW 0IA6RAM (»TTI> 
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FLOW DIAGRAM (OPT!) 
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FLOW DIAGRAM (®PTi) 
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FLOW DIAGRAM (®PTI) 
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FLOW DIAGRAM ( 4>PTI ) 


ENTRY 

$PTC 


INOSKP ! 


PRINT : t, t tl rig, Yg, h, u, v, 

P, q. r, M n , Vq, q* X g , Y g , Z q 


| - -- {^hg>AMAXH J 
,, , YES 

INOARC * F INOARC -T 


I CALL STGTST 


RETURN 


S >0 >' 


INOAPC * 0 

___«■' 

YES ^_ 


PRINT: a, 0 







INDGSR * 0 

- _' 

YES ^_ 

INDFPA = 0 
YES ^_ 

INOFPR * 0^ 

YES L_ 


PRINT : V n 


PR I 


NT . y, <r 


PRINT y, cr 


r 

I 


INOSTE * 0 



PRINT : 6p, ^p,<#.p | 

-T- 

1 



_3 

WRI 

'ST4>P- 


TE : 
HGCTF" 


PRINT : 
Aip, Ayp, Agp 


-(lNOACM«O y - 2 PRINT; a x ,a y ,a x 


:v 


RETURN 


1 INDVYVT • nV_ N ® 

PRiPOT ■ W*r 




PRINT TIME HISTORY 
FR0M LGEAR.VPCS 
TFFS, SACS 


)RY 

9 PRINT F 0c f 

F c i, F ci 


INORMC * 0^ 
I YES 


PRINT; 9, 


\ 
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FLOW DIAGRAM {<J)PTi > 
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(2) CALL LGEAR2 

Initialization after data read in. 

(3) CALL LGEAR3 

Landing gear dynamics computation, Part I 

(4) CALL LGEAR4 

Initial Print. (None in this subroutine). 

29. LGEA3C - LANDING GEAR CALCULATIONS, PART II 

a. Purpose - Continuation of the calculations of the effects of 
ground reaction and landing gear dynamics on the motion of a landing 
vehicle. 

b. Usage - Linkage to LGEA3C is provided by the following statement: 
CALL LGEA3C 


30- SDFLGP - PRINTING OF SDF AND LG VARIABLES 


a. Purpose - The SDFLGP subroutine saves data on tape for plotting, 


prints output, and updates variables for the integration routine. 

b. Usage - Linkage to SDFLGP is provided by the following statements. 
(1) CALL SDFLGP 


Saves data on tape for plotting. 


Input cards required: 


Coiumn: 


1 


o 

L 


IPLT 1 
I SDF 1 
ISTPL1 1 
ISTPL2 1 
IS1PL3 1 
ISTPL4 1 
ISTPL5 1 
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IPLT = 1 denotes that data will be saved on tape unit 13 for plotting. 
ISDF=! denotes that rigid body data will be saved on tape.ISTPL denotes 
data for landing gear k will be saved on tape. 

(2) CALL LGEA6P 

Writes code on tape to designate last record for a particular 
case and also writes EOF. 

(3) CALL LGEAR6 

Prints output of the Landing Gear calculations 

(4) CALL LGEAR7 

Updates landing gear variables that are integrated by the 
integration routine. 

(5) CALL LGEAR5 

Prints titles. (Not used, only a return statement.) 









AFFDI, TR-7M55 
PART IV 


FLOW Qj AG RAM (LGEARI) 
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FLOW DIAGRAM (LGEARl) 



UO 






















FLOW DIAGRAM (LGEAR1) 
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FLOW DIAGRAM (LCEAR1) 
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FLOW DIAGRAM (L GEAR I ) 



NSTRUT NSTRUT NSTRUT 

f TRX 1 Z F TRXj , f TRY ■ 1 ' F TRYj • f TRZ = Z F TR2 

i = I i = I i = I 


f tra 


"Lll R LI2 r L!3 


f TRX 

f trb 

- 

r L2I r L22 r LE3 

• 

f try 

f trc 


r L3I r L32 r L33 


f trz 


NSTRUT NSTRUT NSTRUT 

“TRX * Z m TRM; . “TRY 1 Z “TRY, . “TR2 = Z “TRZ: 

i = I i = l i=l 


M.x 


r lii 

r LI2 

r LI3 


“TRX 

Mty 


r L2I 

r L22 

r L23 

■ 

“try 

Mth ^ 


r L3I 

r L32 

r L33 


“trh 

L J 


NSTRUT 

f XM = F TRA + Z ">j *i o 1 31 . ^YM s r TRB 
i -1 


NSTRUT 

f ZM 5 f TRC + Z "i »'■ °i 33 
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FLOW DIAGRAM { LGEAR I) 
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FL OW DIAGRAM LGEA3C 
(i = 1,2,_ NSTRUT IN Tills SUBROUTINE ) 



RAX j 

*1 



m, 

"i 


Rix 

+ Olj, 

(r Fj 

- Si ) 

RAYj 


8 

/« 


n 2 





Riy 

RAZj 

- 


4 

m. 

"» 


Ri* 

+ °'S3 

(rFj 

- Si 1 

^oi 

= 

( Xg - 

Rg R + 

> 

X 

RG S + (Zg + RAZj ) 

RG,j 



**, 

- 

(Xg - Rg R + 

RAXj) R6|| + (Zg + RAZj 

)R6 

13 



Y Rj 

r 

Yg + 

RAYj 









FOR i = 1,2,-, NSTRUT 
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FLOW DIAGRAM LGEA 3 C 



Rt >*i 5 “ Si a| 33 " (rp i “ q °l 3 | + PR| y “ qRi x 


R DXGj 


r • 1 

xg 


'l 

m, n, 


R DX j 

R DYGj 

s 

Yg 

+ 

*2 

* n 2 n 2 

S 

R D Y j 

R DZG j 


Zg 


'3 

m 3 n 3 


R DZj 

OH 


VTX j = RG,| RDXGj + RG, 3 RDZGj -t- WTj RI 2 2i < ro j — Sj ) 
VTYj » ROYGj - WT } RI 2 , j ( ro ; - 8;) 

VTZj * RG 3 j RDXGj + RG 33 RDZGj 
VGPTj « ( V*TXj -*• V*TYj ) i 
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FLOW DIAGRAM (SDFLGP) 


SOFuGP ' 
ENTRY , 


IPLT = 0 
I NO 


(RETURN 


SAVE THE FOLLOWING 
VARIABLES IN ARRAYS: 
p T2.Sf2> S 2> p 2. p 22. 

s 2 ,s 2 ,s 2 , s 22 ,s 22 , 

s 22: M A2> **^T2• w T2 


LASTPT 3 0 
IPPT *IPPT +I 
(SAVE t IN 
AN ARRAY) 


ISTPL3 * O y -— 
NO 

__jr_ 

SAVE THE FOLLOWING 
VARIABLES IN ARRAYS: 

r T3. Sf3.§3. p 3. p 23, 
S3, S3,S3, S 2 3 . S23, 
S 23> m A3> “T3- w T3 


SAVE THE FOLLOWING 
VARIABLES IN ARRAYS: 
L m, M m, N m. 1 , ®p, 

'f'p* ^*P' p *m> a z« 

Kg , Yg, a„, Xg 


ISTPLI* 0 

~ ~™Tno 


^ ISTPL< • O J— 

nio 

_ ][ _ 

SAVE THE FOLLOWING 
VARIABLES IN ARRAYS: 

p T4> S F4>&4> p 4» p 24> 
S 4 , S 4 , S 4 , S 2 < j, S 24 , 

S Z 4 ,M A 4 ,u; T4 , w t4 


SAVE THE FOLLOWING 
VARIABLES IN ARRAYS: 

P TI > Sp|, 8 |, p |. p 2 (. 
S|» S|, S„ S 2( , S 2 |, 

s 2|. m AI> ***TI. W TI 


ISTPL2* O 
YES 


^ ISTPL5 * 0 J - 

“'"no 

'' 

SAVE THE FOLLOWING 
VARIABLES ON ARRAYS: 

P T5> S F5- ®5« p 5* P, 5, 

s 5 ,s 5 ,s 5 , S 26 ,s 25 , 

®25• M A5» ‘‘'TS- w T5 
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FLOW DIAGRAM ISDFLfiPy 
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FLOW DIAGRAM (SOFLGP) 


WRITE ARRAYS OF THE 
FOLLOWING VARIABLES 
CN TAPE <IPPT POINTS 
FOR EACH VARIABLE I: 
f T 5. Sp5, 85 ,Pa.Pjj, 
S5,Sg,S 5 ,S25, ^ 25 " 
S 25 , M A5 , w r5 , w T5 


IPPT'O 


/TIME LEFT 
ON COMPUTER 
>20 SECONOS ? 



YES 


NO 


LASTPT» I 

WRITE LASTPT ON TAPE 
LTP = 0 

WRITE END-OF-PILE 
CN TAPE UNIT 13 




LTPT 1 0 
LASTPT = I 

BACKSPACE TAPE UNIT 
13 ONE RECORD 


WRITE LASTPT ON TAPE 


rV/— 

IRETURNH*-LASTPT * 


WRITE END-OF-FILE 
ON TAPE UNIT 13 
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FLOW DIAGRAM ( SDR LG P) 


ENTRY\ 

IlgearJ 


INOLG•0 


WRITE :"LGEAR“ 


INDLG * 2 


WRITE NAMES AND 
VALUES OF THE 
FOLLOWING VARIABLES: 
s i. Pi. P2i. PtI.Srj, 

[ s Fi i Oj. Fczi. /*i. 
V GPTi> P'i'Rxi. F TRyi. 
[TRlj. M Ai , M 0 j, 

88), S,, Sj, S,, 

^2i. S2j, S 2j , w Tj , 
w Ti < • ' I. .. NSTRUT) 


write names and 

VALUES OF THE 
FOLLOWING VARIABLES: 

F Ti. r*i. F TR,i, 
F TRyi - F T f}|r j, 88 j, S j 
NSTRUT) 


WRITE NAMES AND i 

VALUES OF THE ] 

FOLLOWING VARIABLES: 
f tra. f T rb. f T rc.m t », 

F,,,,. F ym , 

F im* L m> 


{RETURN 


RETURN) 


RETURN) 


f ENTRY \ 
Ilgear7; 


XNDLG'O 


INDLG »2 


UPDATE THE FOLLOWING 
VARIABLES FOR THE 
INTEGRATION ROUTINE: | 
^2•. Szi 

• ' . .. NSTRUT) 


UPDATE THE FOLLOWING 
VARIABLES FOR THE 
INTEGRATION ROUTINE: 

^ i. Sj, u) j-j 

■ I, .... NSTRUT) 


RETURN 


{RETURN 
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SECTION V 
FORTRAN EXTENDED 
OVERLAY (1,0) 


1. T0LAN1 - MAIN PROGRAM FOR OVERLAY 

a. Purpose - A main program is required by Fortran Extended for 

each Overlay. 

b. Usage - Linkage to Overlay (1, 0) is provided by the following 
statement: 

CALL OVERLAY (T0LA1, 1,0) 

c. Flow Chart 


call diroda 
call DiAiDA 
CALL DIR2DA 


CALL 

READ 



CA 

TFf 

LL 

>S2 

1_1 
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2. TFFS2 - SET UP TABLE ROUTINE 

a. Purpose - To call TSRCH to set up table for various subroutines. 

b. Method - When this routine is called by a certain subroutine, 

it will in turn call TSRCH in order to search the directory for the exact 
BCD name required by the particular subroutine that called this routine. 

c. Usage - Entry is made to the routine with the following statements 

(1) CALL TFFS2 

Sets up tables (or provides the table subscripts) for the 
subroutine TFFS. 

(2) CALL SACS2 

Sets up tables (or provides the table subscripts) for the 
subroutine SACS. 

(3) CALL VPCTAB 

Sets up tables (or provides the table subscripts) for 
the subroutine VPCS2. 

(4) CALL OPTTAB 

Sets up tables (or provides the table subscripts) for 
the subroutine SDF2 (0PT2). 

(5) CALL LGTAB 

Sets up tables (or provides the table subscripts) for 
the subroutine LGEAR. 


137 



AFFDL-TR-71-155 
PART IV 

3. READ - INPUT ROUTINE 

a. Purpose - To provide a general method of reading a variable 
field data card and assigning variable length table. Data may be read 
into symbolic locations in memory. 

b. Method - Decimal, Octal, and Integer numbers are converted to 
binary integers. BCD information is stored in six-character words. 

c. Card Formats - The variable name punched in Columns 1-6 is the 
location into which the first data word is loaded. The variable field 
information is for relocation. A fixed point integer punched anywhere 
in the field (67-72) will be treated as a subscript to the variable name 
punched in Column 1-5. Negative integers punched in columns 67-72 will 

be in violation of the subroutine. The first blank character found in the 
card to the right of Column 12 terminates loading from che card. One 
exception to this is BCD data. Nine 6-character words may be loaded 
including blanks. 

The general character of the data tc be loaded is determined by a 
three letter pseudo-operation punched in columns 8-10. The pseudo¬ 
operations are: DEC or blank, 0CT, INT, BCD, and TRA. The pseudo¬ 
operation TRA is a method of exit from the subroutine. 

d. Decimal Data - Decimal data beginning in column 12 and ending in 
column 66 is converted to binary and loaded into the symbolic location 
punched in column 1-6 subscripted by the integer punched in column 67-72. 
Signs are indicated by + and - preceding the number. All unsigned 
numbers are treated as positive. If either the character E or . or 
both appear in the decimal data word, the word is converted to a floating 
binary number. The decimal exponent used in the conversion is the 
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number which follows immediately after the character E. This number may 
have a + or - sign preceding it. If the character E does not appear 
the exponent is assumed to be 2 ero. If a decimal point does not appear 
it is assumed to be at the right of the numbers, unless it is the only 
word or the first word on a card; then it is assumed to be an integer. 

All the examples below are equivalent. 

(1) 12.345E03 

(2) 12.345E+03 

(3) 12.345E3 

(4) 12345E00 

(5) 12345. 

(6) 1.2345E4 

(7) 1234500E-02 

(8) +1234500E-2 

Note that in the examples above all decimal words have decimal points. 

If the frst word c.n a card, or if it happens to be the only word on a 
card, does not contain a decimal point, the word will be converted to 
binary integer. 

e. Octal Data-0ct - The Octal data is loaded the same as decimal 

data but must have 0CT punched in columns 8, 9, and 10. All data is converted 
to binary with binary point assumed at the right end of each word. 

f. Hollerith Data - BCD - Hollerith information is loaded from 
columns 3 through 66 and assigned consecutive locations for every 6 
characters. A maximum of nine 6-character words may be punched on 
an.v one card and the number of words must be punched in column 12. 

A subscript may also be punched in columns 67-72. 
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g. Transfer - TRA-- The purpose of the TRA card is to transfer 
control from the subroutine back to the main program. TRA must be punched 
in columns 8, 9, and 10. The subscript field is not used. A REWIND may 
be punched beginning in column 12. Only the R is checked and the only 
use is for the rewind of a data tape. 

h. Integer - INT - Integer data begins in column 12 and ends in 
column 66. INT is punched in column 8, 9, and 10. It may be relocated 
with respect to the BCD name by punching a subscript integer in 
columns 67 through 72. If only one data word is punched per card, 
columns 8, 9, and 10 may be left blank. 

i. Error Messages - A message is written on the output tape 
describing the type of error encountered. If an error is encountered, 
execution of the case is deleted and the subroutine only searches for 
other possible errors in the data. The following error messages are 
possible. 

(1) Symbol not in directory. 

(2) Column 12 is blank. 

If a bad pseudo-operation is punched in columns 8, 9, and 10 the 
subroutine will treat it as decimal data. 

All checking for redundancies, end of tape, format errors, etc., 
ns handled by FORTRAN system input/output routines. 

j. Usage - Linkage to the routine is made by the following statement: 
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k. Data Preparation - The first card expected by READ is a 
STCASE TAB with the S beginning in Column 1 and TAB punched in 8, 9, and 10. 
Following this card is a set of cards which define the table sizes necessary 
for that case. 

Example: CTAB01 10 

CTAB02 20 

On the above example, CTAB01 is punched beginning in Column 1. The 
numbers 10 and 20 are punched in Column 12 and indicate the number of 
machine cells necessary for that table. Any number of tables may be 
assigned as long as the total number of machine cells does not exceed 600. 
Follow all table assignments with a TRA punched in Columns 8, 9, and 10. 

The next data required by the subroutine is a STCASE and followed by 
any combination of 0CT, BCD, INT, and TRA cards. 

4. DSERCH - DIRECTORY SEARCH ROUTINE r 0R SUBSCRIPTS 

a. Purpose - To provide a method of searching the directory to 
find the subscript corresponding to a BCD argument. 

b. Method - The routine searches the directory for the exact BCD 
name required by the argument. When an equal compare has been found, 
the corresponding subscript is returned as a fixed point integer. 

c. Usage - Entry is made to the routine with the following 
statement: 

CALL DSERCH (SYM, L0C, IER) 

where 

SYM = BCD name being search for. 

L0C = The location in core 
IER = Error Code: 

IER is set to 0 if the BCD argument does compare and 
is not a table name. 
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IER is set. to -1 if the BCD argument does compare and 
is a table name. 

IER is set to +1 if the BCD argument does not compare. 

5. TABRE - TABLE DIMENSION SUBSCRIPT ROUTINE 

a. Purpose - To compute subscripts such that the table dimension 
requirements may be variable. 

b. Method - Uses input data prepared by the user to compute 
subscripts for variable table assignments. 

c. Usage - Entry is made via the statement 

CALL TABRE 

When the READ subroutine processes the control card: STCASE TAB 
then the subroutine TABRE is called by the READ subroutine. STCASE begins 
in Column 1 and TAB punched in Columns 8, 9, and 10. Following this card 
will be the cards requesting table sizes, which are read by the TABRE 
subroutine. As an example the following cards may be read: 

TTAB10 10 

ATAB01 2 

TTAB10 and ATAB01 punched beginning in Column 1 and the required machine 
cells (10 and 2 in this case) punched beginning in Column 12. Anything 
punched past Column 15 will not be used. After all table assignments, 
a TRA should be punched in Columns 8, 9, and 10. The following error 
messages may be printed: 

1. " Symbo l does not exist in table list." 

2. "Total table sizes exceed N, change maximum total table sizes to NN 
in subroutines TLU, HIH0, TFFS, AER0, and AUXR2," where NN is the 
required and N is the maximum. 
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6. TSRCH - TABLE SUBSCRIPT SEARCH ROUTINE 

a Purpose - To provide a method of searching the directory for 
table subscripts. 

b. Method - The routine searches the directory for the exact BCD 

name required by the argument. When an equal has been found the corresponding 
subscript is returned as a fixed point integer. 

c. Usage - Entry is made to the routine with the following statement: 
CALL TSRCH (SYM2, L0C2, N2, IER) 

where 

SYM2 = BCD name of argument 
L0C2 = Location of first subscript 
N2 = NuPiber of sequential subscripts to return with 
IEK = Error Code: 

IER is set to +1 if the BCD argument does not compare 
IER is set to -1 if the BCD argument does compare. 
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flow diAprak 1 ;tabre.) 
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7. ROUTINES CALLED BY READ SUBROUTINES 

The following are subroutines called by the READ subroutine for which 


flow charts are not provided. 
NAME OF SUBROUTINES 
DIPLAC 
READA 
STORE 
WRCARD 
PACKRR 
RITE 


LINKAGE 

CALL DIPLAC (RA1, INC, BLANK) 
CALL PEADA (IBCRW) 

CALL'STORE (N, INC, INX, STAPE) 
CALL WRCARD (MSG) 

CALL PACKRR (II, 12, NNN) 

CALL RITE (IFI, FO, JJ) 


8. DlftOQA - INPUT DIRECTORY, PART I 

A routine of all the variables in data statements that may be read 
by the input routine and their corresponding subscript lccation in 
COMMON/DIRCOM. 

9. DIRIDA - INPUT DIRECTORY, PART II 
Continuation of DIRO. 


10. DIR2DA - INPUT DIRECTORY, PART III 
Continuation of DIR1. 
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SECTION VI 
FORTRAN EXTENDED 
OVERLAY (2, 0) 

1. T0LAN2 - MAIN PROGRAM FOR OVERLAY 

a. Purpose - A main program ■'s required by Fortran Extended for 
eacn Overlay. 

b. Usage - Linkage to overlay (2, 0) is provided by the following 
statement: 

CALL OVERLAY (TOLA, 2, 0) 

c. Flow Chart 
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2. AUTS - AUTOPILOT AND CONTROL SYSTEM 

a. Purpose - The AUTS subprogram calculates the control responses 
that are needed by the aerodynamic forces and moments subprogram. 

b. Usage - Linkage to AUTS is provided by the following statements: 
CALL AUTS - Compute the control responses 

CALL AUTSPR - Compute the control response and print information 
as indicated by input parameters. 

3. FLARE1 - AUTOPILOT FLARE 

a. Purpose - The FLARE subroutine calculates the actual selected 
touchdown conditions, time to touchdown for X and Y constraints, required 
accelerations to meet touchdown conditions, airspeed flight path angle 
relative tu runway, desired thrust, and angle of attack. 

b. Usage - Linkage to FLARE is provided by the following statements: 
CALL FLARE1 (IPR) Computes all that is mentioned above. 

IPR is the print indicator. If IPR = 0, 
the print routine is called, 

CALL FLARE2 (IPR) Only computes airspeed flight oath angle 

relative to runway, desired Euler roll angle 
and desired angle of attack. 
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FLOW DIAGRAM (AuTS) 
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FLOW DIAGRAM (AUTS) 



O' 



IAP « 1 

v,*v g -v d 

1 


'|v,j : 

1 

>* 

< 

A 

1 

YES 

r_ 


hi© 


v Od * lV d CO* «gg - Stga»i* 

■M Yg, I* + (V d tin « 6S -Zg* >*)’'* 
XjJ * TAN - I [(- Vj tin - Zg« I 

/ V d CO* «qs ~ Xg» I] 

Qr * sp^ca* 

C LR * ^ g REF / ,0 2 S * 

(THIS IS A FIRST APPROXIMATION 
FOR C LR ) 

SWT2 » . FALSE. 

9<J * O. 

I ’ O 


( OBTAIN 8o W ■ f («<». SbVcs)’ 
CALL SACSIO 

(OBTAIN « d -f(cLR.S 8 V ES ) 
AND c 0R • f(« d n 
CALL SACS8 
(ITERATE 5 TIMES) 


NO/'--\YES S~\ 
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4. AUTPR1- AUTOPILOT PRINT ROUTINE 


a. Purpose - The AUTPR subroutine prints output for the autopilot. 

b. Usage - Linkage to AUTPR is provided by the following statements: 

(1) CALL AUTPR1 Write ! AUTS". 


(2) CALL AUTPR2 

(3) CALL AUTPR3 

(4) CALL PR4 

(5) CALL AUTPR5 

(6) CALL AUTPR6 

(7) CALL AUTPR7 


If AUXICP f 0, WRITE 

"AUXILIARY COMPUTATIONS" 

and the following variables: X R , Y R , Z R , 

X R , Y r , Z r , <|y 

If MANL0G f 0, write "GLIDE SLOPE" and the 

following variables: Ve, V ad , yj, Q R , C L , 

R 

a d’ L R’ D R’ TTd ’ H GS* R R* h ea’ h e* h eT 5 h Pa’ 
h PT* ♦d 

If MANL0G f 0, write "FLARE" , and the following 
variables: t^, Xjq, ^TD* 

a • 

X RF . X rf , h RF , h RF , A xr , A hR , a d , T Jd , <J> fl 
If MANL0G t 0, write “LANDING ROLL" , and the 
following variables: ISS, ILR, ICS, BS, Tj, t r 
if PITCHP f 0, write "PITCH AUTO PILOT, " and the 
following variables: 6 c ,, a, a d , 

ot t » <5 

eT’ qd 

If YAWAUP f 0, write "YAW AUT0PIL0T" , and the 
following variables; S rN> 6, 6, B eJ , 6 rd> U> et 

If R0LLAP f 0, write "ROLL AUT0PIL0T" 
and the following variables: 

^e* ^eT’ 
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(8) CALL AUTPR8 


(9) CALL AUTPR9 


(10) CALL AUTP10 


(11) CALL AUTP11 


If THR0AP f 0, write "THR0TTLE AUT0PIL0T" , and 

the following variables N.., I. , i = 1, 2, 

a a i 

.... IN 

If BRAKAP f 0, write "BRAKE AUT0PIL0T" and the 
following variables: , i = 1, 2, 

NSTRUT. Also if BRAKAP f 0 and IBS = 1, 
write the following brake autopilot variables: 

u tr i» ^Ri* w i£i» 1 = 1 » 2 » NSTRUT 

If C0NTRP t 0, write "C0NTR0L RESP0NSE", 

and the following variables: 6„, 5 , 6 , 

q r p 

6 , & r , 6 , N jt N-, i = 1, 2, .... IN. 

Also, if INDICP f 0, write "INDICAT0RS" 
and the following variables: IAP, IC i , 
i = 1, 2, 3, 4; IB i# i = k, 2, 3, 4, 5. 

If ‘1ANL0G f 0, write "FLARE", and the 
following variables: ,, d> 0 rf , T Jd 


5. THAUTS - AUTOPILOT THROTTLE 

a. Purpose - To compute the desired setting and desired thrust for 
each engine. 

b. Usage - Linkage to the throttle autopilot is made by the following 
statement. 

CALL THAUTS (IPR) 

If IPR = 0, the print routine is called to print the throttle 
autopilot output variables. 
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FLOW DIAGRAM (AUTPRI) 



^MANL0g « 0^- 


YES 


NO 


WRITE.- LANDING 
ROLL “ 

WRITE LANDING R$LL 
VARIABLES: ISS, 
ILR, ICS, IBS, 

T I.*r 






NO 

I WRITE : "PITCH 

1 AUT#PILI>T " 

1 WRITE PITCH 
AUT0PIL0T VARIABLES: 

8qN, <*•, <i 
a *T. 8 qd 

. «*d. 




0 




YAWAUP 


-y« 

P - 0 y ^{return) 

no v^y 


WRITE : “YAW 
AUT$PIL$T" 

WRITE YAW AUT$ - 
PIL^T VARIABLES: 

& r N> 0, &, $*T> 
fer d' *•- i'eT 
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6. ENGREV - AUTOPILOT ENGINE REVERSE LOGIC 

a. Purpose - ENGREV is a function that determines if an engine 
is to be reversed. 

b. Usage - Linkage to ENGREV is provided by the following statement 
Y = ENGREV (IR, N, N g . N,_ R ) 

where 

IR = engine reverse capability 
N = engine throttle setting 
N„ = engine throttle constraint 

D 

^LR = reverse throttle setting for reverse signal on landing 

7. ENGFL - AUTOPILOT ENGINE FAILURE LOGIC 

a. Purpose - The ENGFL subroutine contains the logic to determine 
if an engine has failed or not. 

b. Usage - Linkage to ENGFL is provided by the following statement: 
CALL ENGFL 

8. CTENGL - AUTOPILOT COMMON TWO-ENGINE LOGIC 

a. Purpose - The CTENGL routine calculates the distribution of a 
totei dosireci thrust load, T , between two engines, A and B, as T.. 
and T dg respectively. It al.o determines the corresponding desired 
throttle settings, N dA and N dg , by calling CENGL. twice. 

b. Usage - Linkage to CTENGL is provided by the following statement: 
CALL CTENGL 
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9. CENGL - AUTOPILOT COMMON ENGINE LOGIC 

a. Purpose - The CENGL routine calculates the variable N^ which 

c 

is required by throttle autopilot (THAUTS) 

b. Usuge - Linkage to CENGL is provided by the following statement: 
CALL CENGL 

10. TFFS8 - THROTTLE SETTING SEARCH ROUTINE 

a. Purpose - Search through ranges of throttle settings to find the 
throttle setting TN that corresponds with current thrust TC and Mach 
number. 

b. Usage - CALL TFFS8 (TC, TN) 

11. TFFS9 - COMPUTE THRUST as F (MN, TN) 

a. Purpose - Compute thrust as a function of current Mach number 
and designated throttle setting. 

b. Usage - CALL TFFS9 (THRSET, THRUST) 
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FLOW DIAGRAM t TFFSI) 


GENTRY' 
i TFFS8i 


INDTFF'O 

I YES 

r Tn «o. 


RETURN 


LOCT * LOC(I) -I 
WHERE LOC(I) 5 
LOCATION OF THROTTLE 
SETTING VALUES, (N), 
IN THE THRUST TABLE 
I = 0 


TN' 

■-2. - 


NO 

^ TC> TH2 

, r 

YES 


: 

\ 

THI = 

TH2 


IfcITlOW 
I YES 



TN * C(K -1)+( THI 
-TCI*(C(K) 

-C (K-IH/ITHI 
-TH2) 


I 5 I + I 

K • LOCT + I 
TH2»f (C(KJ,M n ) 
WHERE C IS THE THRUST 
TABLE AND N IS THE 
SUBSCRIPT LOCATION 
OF N(THROTTLE 
SETTING! 



I NO 
TC>TH2 


C(K) > -2 


ES 


j NO 
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SECTION VII 

PLOT TAPE GENERATING PROGRAM (PLTSDF) 

a. Purpose - To generate a plot tape from a data tape that was 
generated by the Takeoff and Landing Analysis Computer Program. The plot 
tape is plotted by the CALCOMP plotter. 

b. Usage - After a plot tape has been generated by the TOLA program 
(see subroutine SDFLGP for data setup), this program is turned in as a 

s parate job. 

c. Input: (1) TAPE - Data generated by TOLA program. 

(2) CARDS - The following may be read in as input using 
the NAMELIST feature of F0RTRAN EXT. If the value of any variable is the 
same as its nominal value, it is not necessary to read it. 


Name Definition Nominal Value 

NCASES = Number of sets of data or cases to be plotted 1 

ISDFR = 1 Rigid body data is stored on tape 1 

= 0 Rigid body data is not stored on tape 
ISDF = 0 Do not plot rigid body data 0 

= 1 Plot rigid body data 

ISTPR(l) = 1 Landing gear data for gear 1 stored on tape 1 


= 0 Landing gear data for gear 1 is not stored on tape 
ISTPL(I) = 0 Do not plot landing gear data for landing gear 1 0 

= 1 Do not plot landing gear data for landing gear 1 
IL = 1 Do not plot lower chamber pressure and 2nd Piston 0 
plots (2nd piston was not uspd) 

= 0 Do plot lower chamber pressure and 2nd piston plots 
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Name 

TFIRST = 
TLAST = 


PLTINT = K 
FCTR = 


XL 

YL 


Definition Nominal Value 

Trajectory time to begin plotting 0. 

Trajectory time to stop plotting ( Note= 0. 

if both TFIRST = TLAST = 0., the entire 
time history on tape will be plotted) 

Plot every point 1 

The current factor all coordinates are 
multiplied by. (That is, the plot is made 
larger or smaller if FCTR is greater than 1. 
or less than 1. For example, if wish plot 
to be 25% of the original size, let FCTR=.25) 

Length of X - Axis in inches 7.2 

Length of Y - Axis in inches 5.0 


Examples of input data cards 

Example No. 1: Plot rigid body variables and landing gear variables 
for gears 1, 3, 5. Plot every point and plot entire time history. Assume 
rigid body and landing gear variables for gears 1, 3, 5 are stored on tape. 
EOR (7/8/9) 


SINPUT ISDF=1,ISTPR=1,0,1,0,1,ISTPL=1,0,1,0,1$ 

E0J ( 6/7/8/9) 

Example No. 2: Plot landing gear variables for gear No. 3. Plot 
every other point from time = 4. to time * 10. seconds. Assume rigid 
body and landing gear variables for gears 1, 2, 3, 4> 5 are stored on tape 
for time = 0. to 20. seconds. 

EOR (7/8/9) 

$INPUT,ISTPL(T)=1,TFIRST=4.,TLAST=10.,PLTINT=2$ 

E0J (6/7/8/9) 
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Example No. 3: Plot rigid body variables and landing gear 
variables for gear No. 5. Plot every point and plot entire time history. 
Assume rigid body and landing gear variables for gears 1, 3, 5 are stored 
on tape, with the size of graphs to be 50% of the original size where the 
original size of the X-Axis is 8 inches and the Y-Axis is 6 inches. 

EOR (7/8/9) 

$ INPUT,ISDF=1,ISTPR-1,0,1,0,1,ISTPL(5)=1,XL=8.,YL=6.,FCTR=.50$ 

EOJ (6/7/S/9) 


d. Output 

TAPE. A plot tape is to be plotted by the CALC0MP plotter. The 
plot tape may be generated to give the following plots: 

(a) Rigid body plots 

Time (sec) Vs Ground Reaction roll moment (lb>^0 
Vs Ground reaction Ditch moment (Ib-ft) 

Vs " " yaw moment (Ib-ft) 

Vs Pitch rate (rad/sec) 

Vs Yaw angle (Deg) 

" Vs Roll angle (Deg) 

Vs Z-Axis ground reaction (lb) 

Vs Z-Axis acceleration (ft/sec^) 

Vs Mass Center sink rate (ft/sec) 

Vs Mass center altitude (ft) 

Vs Cross range (ft) 

Vs Mass Center altitude (ft) 

Vs Cross range (ft) 

Vs X-Axis acceleration (ft/sec^) 

Vs X-Axis Velocity (ft/sec) 


Downrange (ft) 

ii 

Time (sec) 
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(b) Landing gear plots for each gear (1) 

Time (sec) Vs ground reaction along strut for gear No, 1 (lb) 
Vs strut resistance for gear No. 1 (lb) 

" Vs tire deflection for gear No. 1 (ft) 

2 

" Vs upper chamber pressure for gear No. 1 (lb/ft. ) 

" Vs lower chamber pressure for gear No. 1 (lb/ft.^) 

" Vs strut acceleration for gear No. 1 (ft/sec ) 

" Vs strut velocity for gear No. 1 (ft/sec) 

" Vs strut displacement for gear No. 1 (ft) 

" Vs 2nd piston acceleration for gear No. 1 (ft/sec ) 

" Vs 2nd piston velocity for gear No. 1 (ft/sec) 

" Vs 2nd piston displacement for gear No. 1 (ft) 

Vs wheel axle moment for gear No. 1 (lb/ft) 

2 

Vs tire angular acceleration for gear 1 (rad/sr„ ) 

" Vs tire angular rate for gear No. 1 (rad/spc) 

(c) The following routines are calle*. by SDFPLi: 


PL0T 

SCALE 

PL0TE 

PL0TS 

LINE 

READ 

FACTOR 

SYMB0L 

BKSFIL 

AXIS 




(d) The appropriate control cards are required after the job card 

to mount tape xxxx from TOLA on tape unit 13 and to mount tape yyyy for 

FLTSDF output (556 BPI) on tape unit 7. 

(e) If TOLA and PLTSDF are run on the same job, it ’ c not necessary 

to use tape unit 13 for storing data. In this case disk storage may be 

used. 
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